Typeset with jpsj2.cls <ver.1 .2> 



Letter 



Evidence for Octupole Order in Ceo jLao 365 from Resonant X-ray Scattering 

Hiroaki Kusunose * and Yoshio Kuramoto 
Department of Physics, Tohoku University, Sendai 980-8578 

The azimuthal angle dependence observed in the resonant X-ray scattering in phase IV of Ceo yLao sBg 
is analyzed theoretically. It is shown that the peculiar angle dependence observed in the E2 channel is 
consistent with the r5„-type octupole order with principal axis along (ill) and equivalent directions. Under 
the assumption that the four equivalent octupole domains are nearly equally populated in the sample, 
the observed angle dependences are reproduced by calculation for both era-' and crn' polarizations. The 
calculation for various symmetries of order parameters excludes unambiguously other order parameters 
than the r5„-type octupole. 

KEYWORDS: octupole moment, resonant X-ray scattering, rare-earth hexaboride 



Octupole orders have attracted growing interest recently in 
relation to unusual properties found in NpOa,'"^ URuaSia,'* 
and Ce^Lai-iBe.^"^ In particular, Ce^^Lai-j^Be exhibits a rich 
variety of the H-T phase diagram, called I-IV.^ '' The impor- 
tance of the octupole degrees of freedom was first recognized 
by Sakai and co-workers'"'" in resolving the inconsistency 
between NMR and the neutron scattering measurements in 
the antiferro-quadrupole phase II. Soon after their work, pos- 
sible relevance of the Fsu-type octupole to the magnetic phase 
III and III' was pointed out.^ For consistent understanding of 
the overall phase diagram, this type of octupole was also sug- 
gested as a hidden order parameter of phase IV. The predic- 
tion of the accompanied lattice distortion along the [111] axis 
is indeed observed.'^ The model of the Fsu-type octupole or- 
der explains most of the characteristic behaviors in phase IV, ^ 
such as the cusp-like behavior in the uniform susceptibility,'' 
no significant magnetic reflection in neutron scattering,'^ and 
the huge softening of C44 mode in the strain susceptibility'"* as 
well as the spontaneous lattice distortion along [111].'- How- 
ever, there has been no direct evidence from diffraction exper- 
iment. 

Resonant X-ray scattering (RXS) measurement has been 
performed recently to obtain a direct evidence of the order 
parameter symmetry in phase IV. '^ The RXS has already 
proven its efficiency in various systems such as LaMnOa,'^ 
DyBaCa'*"^ and NpOa.^'^^ The measurement was done for 
Ceo.7Lao.3B6 below Tiv = 1.5 K, using the L2 absorption 
edge in the Ce atom.'^ They measured the azimuthal an- 
gle dependences of the superlattice Bragg reflections at ^ = 
(3/2,3/2,3/2) with the rotation axis normal to the [111] 
surface. Both the non-rotated (cr-cr'), and the rotated (cr- 
n') polarization channels have been measured at the electric 
quadrupole (E2) resonance. In the cr-cr' channel, the scatter- 
ing intensity exhibits six-fold oscillation, which indicates the 
occurrence of the electronic order with three-fold symmetry 
along [111]. In contrast, the cr-n' channel shows three-fold 
oscillation, which remains to be clarified. It is not clear to us 
how the fitting in ref. 19 is justfied microscopically. 

In this Letter, we demonstrate that the F5„-type octupole 
order can reproduce the azimuthal angle dependences both in 
the cr-cr' and cr-n' channels. The three-fold oscillation is ex- 
plained in terms of domains whose principal axis is other than 



[111]. These domains give the scattering amplitude which has 
both three-fold oscillatory component and the axially sym- 
metric one. As a result, their cross product contributing the 
intensity gives the three-fold symmetry. We adopt a concise 
formalism for the RXS amplitude proposed by Lovesey and 
Balcar,^"'^' which assumes the spherical symmetry in the in- 
termediate states. Our results of calculation are compiled into 
tables, which allow us to analyze various multipole orders 
quite easily. With an appropriate rotation to obtain domains 
with different principal axes, we summarize the azimuthal an- 
gle dependences for all possible order parameters. The result 
excludes unambiguously the other symmetries than the F5,, 
octupole. 

The RXS amplitude per unit cell is given by the standard 
formula:^"'^' 

m 

where we have neglected the energy distribution of the in- 
termediate states (m), and of their width F. In the RXS, the 
energy of the incident photon fio) is tuned close to the absorp- 
tion edge A of the relevant atom. The intensity is proportional 
to 

^{f\e' ■P\m){m\e-P\i) 

+ {f\Tv{X' ■Q)\m){m\Tx{X-Q\i), (2) 

where the first term in the right-hand side comes from the 
electric dipole (El) transition, and the second term from the 
quadrupole (E2) transition. Here we have defined the dipole 
and the quadrupole operators in an atom, P" - e and 
Qap - s '"n'^/2> and the matrix made from the photon 
wave and polarization vectors, Xap - k"^ jl. The prime rep- 
resents quantities concerning the scattered radiation. 

Following Lovesey and Balcar,^"'^' we assume that the in- 
termediate states can be described by an angular momentum 
of a core hole. Then, the summation in eq. (1) can be done 
with the help of the Wigner-Eckart theorem, and we are left 
with two scalar products of the spherical tensors; 

X<=4'.zf, 
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Y,T^\M,M')K't{e',e), (3a) 



q=-p 



Zfl = Z^f Z TT{M,M')Hf\Q',Q), 

p=0 q=-p 



(3b) 



where T*''^ represents the multipole degrees of freedom in 
the / shell. It is related to the spherical harmonics as T^^ = 
^jAnlilp + Vjr^Ypq where r is the radial coordinate. The oper- 
ator equivalents of T^\m, M') can be obtained by replacing 
the coordinate r with its symmetric product of the total an- 
gular momentum /, e.g., xy^ — > {JxJ^. + JyJxJy + JyJx)/^- 
The information of k and e is summarized concisely in the 
corresponding spherical tensors, 
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Fig. 1. The coordinate system corresponding to the RXS experiment. The 
azimuthal angle dependence is obtained by rotating the r = (x, y, z) co- 
ordinate system around the z axis relative to the photon k vector, i.e., k is 
in the yz plane att// = 0. The u and n polarization vectors are also shown 
(e^xe„ = fc/|fc|). 



H^g'^ = F Yj Khfi{2a,2/3\pq), 



(4a) 



(4b) 



Table I. The angle dependences of the RXS amplitude in the non-rotated 
(o--cr') and the rotated (cr-n') polarization channels for the El transition. 
The common coordinate systems for the scattering and the crystal is used. 
The definition of the cartesian spherical tensor K^]sq is given in the text. 



ha= Y ^^^(^7' iS\2a), 



(4c) 



7(5=- 1 



where (jimi, j2m2\jm) is the Clebsch-Grodan coefficient. An 
ordinary vector v in the cartesian coordinate is related to the 
1st rank spherical tensor as v'j'j^ = +(v.i + /v,)/ V2 and ■* = 
Vj. For all spherical tensors in this paper, we use the phase 
convention [5^^]* = {-lys'-f^^^ 

The operator / is the axial vector and has an odd parity un- 
der time reversal.^^ As a consequence, within the same man- 
ifold of the angular momentum, an odd-rank tensor opera- 
tor describes a magnetic multipole, which can be coupled to 
an axial tensor with time-reversal odd, such as ie' x e, while 
an even-rank tensor operator represents an electric multipole 
which appears only with the polar tensor with time-reversal 
even. When a spontaneous order of Tg'^ multipole occurs, 
i.e., (Tg''*} + 0, the azimuthal angle dependence of the cor- 
responding will appear in the RXS ampUtude. The scalar 
form of eq. (1) guarantees that arbitrary linear combination of 
T^q ^ couples with the complex conjugate of the same linear 
combination of The scalar property is a great advantage 
in considering the case where the ordered moment is directed 
to a direction different from the z axis. This situation will be 
discussed in more detail later. Since the explicit form of the 
coefficient c^" (n - 1,2) depends on the representation of 
multiplet, i.e., the Russell-Saunders or the j- j coupling, we 
regard it as a fitting parameter. For further details, we refer to 
the papers by Lovesey and Balcar.^'''^' 

The experimental geometry of the scattering is shown in 
Fig. 1 with the coordinate system r = {x,y,z). The azimuthal 
angle i^/ dependence is obtained by rotating the r coordinate 
system around the z axis relative to the photon k vector. The 
k vector is then given by = A:(cos 9 sin i/r, cos 9 cos ijf, sin 6). 
Tables 1 and 11 summarize the Bragg and the azimuthal angle 
dependences of Kq''^ and Hq'^ for the case where the crystal 
coordinate R = (X, Y, Z) is the same as the scattering coordi- 
nate r. Here we have used the real cartesian tensors for ^ > 



cl 



1 



cos 8 sin if/ 
COS 8 cos ifi 



K' 



c2 
(2) 
cl 

s2 



--^cos(2i/f) 



sin(2^) 



1 



sinf 



- sin 8 sin(2(f') 

cos 6 cos if) 
--^ sin(?cos(2(<f) 
--^ cos0sin(;r 



defined as 



V2 
(-1)^ 
y/2i 



(5a) 
(5b) 



A similar definition is applied to H^''^ as well. Due to the 
choice of 60- || x for 1^ = 0, the tfr dependence in the cr-cr' 
channel is the same as that of the spherical harmonics. On the 
other hand, the change of the photon helicity in the cr-n' chan- 
nel gives an additional phase factor n/2 in the i// dependence. 
For an occurrence of spontaneous multipole order, we obtain 
immediately the i// dependence of the RXS amplitude from 
Tables I and II for the El and the E2 transitions, respectively, 
as long as the principal axis of the order parameter is parallel 
to the z axis. 

Now, let us discuss the azimuthal angle dependence of 
phase IV in CeojLaojBf,. The El scattering shows a short- 
range magnetic order, which may be induced from the oc- 
tupole order together with some crystalhne disorder. In the 
following, we restrict our discussion to the E2 transition, since 
it provides the direct information of the octupoles. In order to 
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Table n. The angle dependences of the RXS amplitude for the E2 transi- 
tion. 
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take the [1 1 1] direction as the rotation axis z, we set 



X ^ 
Y 



f 1/V6 1/V3 

-1/V2 1/V6 1/V3 
-2/V6 j 



y 



(6) 



Then, [112] direction is parallel to y, and k at {// = is in 
the yz plane. The origin of i/f is the same as the one used by 
Mannix et al. 

The most promising candidate of the order parameter for 
phase IV is the Fs^-type octupole with the [111] principal 
axis: 



Vl5 



+ Z(X^ - Y^) 
In the scattering coordinate, it has the form 



^2-{X-Y){Y-Z){Z-X). (7) 



(8) 



Thus, from the row of //^j' in Table H, we obtain immediately 
the azimuthal angle dependences of the T^^^^^ order as 



l4«'"(^^)l' = Y^sin2(2e)cos2(3^ir), 



(9a) 



Table HI. The angle dependence of the RXS intensity for possible multi- 
poles in three-fold axis. The rotation axis is [111], and k is in the [111]- 
[112] plane ati/r = 0. 
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l/5„ (^^')r = [cos 9 + 3 cos(30)]^ sin^(3^). (9b) 

Similarly, we derive the angle dependences of other multi- 
poles with [111] as the principal axis. The expressions of the 
multipoles with components p-T are given by 



7^11111 
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n(2)/ 



2- 5g 

p[lll] 

3- 2« 

rllll] 
' 3-4m 

p[lll] 

-4-4g 



^rf3>(.).^rf(r). 



xliii] 

4-5g 



' c3 
2VI5 



V2T 



(10a) 
(10b) 
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9 u ^ ' 9 

and their angle dependences are summarized in Table III. 

The averaged scattering intensity from the other domains is 
calculated as 

1 

144 

\f*r\-n't= ' 



\fs^'''\(rcr')f - sin\20[cos\3f) + 2], 



(11a) 



4gQg-[5(cos(30-5cosr 
- (3 cos(36') + cos 6) cos(6(fr) 
-I- 32 '\^(sin(36l) - 7 sin ff) cos^ 6'cos(3i/r). 



(lib) 



Among all multipoles accessible by the RXS, only the r5„- 
type octupole yields the maximum at i// = in the cr-cr' 
channel, which is consistent with the observed oscillation'^ 
as shown in Fig. 2. In the cr-n' channel, on the other hand, the 
contribution from the other three equivalent domains is much 
larger than that from [111] domain, and the threefold compo- 
nent appears with the maximum at = 0. To the contrary, 
the [111] domain predominates over the other domains in the 
cr-cr' channel. 

Finally, we show the result of the fitting for the observed 

data taken from IVIannix et al}'^ in Fig. 2. We have used the 
fitting formula for the scattering channel a (- cr, n): 

la = A„ [w\fll''\at + (1 - w)|/5f' "(a)|2] , (12) 
with the scattering angle = -39° from the incident pho- 
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Fig. 2. The azimuthal angle dependence. The symbols are taken from Man- 
nix et at The dotted and the dashed lines show the intensities from the 
Tsu-type octupole order described by eq. (12) in the text. 



ton energy and the scattering wave vector ^.'^ We assume 
that the four equivalent domains are equally populated, i.e., 
w = 1/4. The calculated intensity is insensitive to slight 
change of w. The intensity factors are chosen as A,r,T' = 120, 
Acrjii = 70. If the ideal experimental setup were achieved, the 
intensity factors in the two scattering channels should be the 
same. The discrepancy of the two intensity factors in the fit- 
ting may be ascribed to slightly different setup for different 
polarizations, and/or extrinsic background from non-resonant 
contribution. The fitting with these intensity factors repro- 
duces semi-quantitatively the observed data. Thus, we con- 
clude that the r5„-type octupole in three-fold axis with the 
four equivalent domains is the order parameter of phase IV in 
Ceo.vLao.sBg. 

In sunnmary, on the basis of the formalism by Lovesey and 
Balcar, we have derived azimuthal angle dependences for all 
multipole orders. The appropriate rotation, which transforms 
original expressions of irreducible tensors into the scattering 
coordinate system, allows us to extract the azimuthal angle 
dependences for all possible candidates of the order param- 
eters for phase IV in Ceo.7Lao.3B6. In the non-rotated cr-cr' 
scattering channel, only the r5„-type octupole in three-fold 
principal axis yields consistent six-fold oscillatory behavior 
with the maximum at = 0. On the other hand, in the rotated 
cr-n' channel, the contribution from the other three equivalent 
domains predominates over the contribution from the [111] 
domain, showing three-fold oscillation with the maximum at 



3) 



4) 
5) 
6) 
7) 



t// = 0. Therefore, we conclude that the Fsu-type octupole in 
three-fold axis with the four equivalent domains is the order 
parameter of phase IV in Ceo.7Lao.3B6. 
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